Infrared reflectivity of the phonon spectra in multiferroic TbMnOg 
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We measured the temperature dependent infrared reflectivity spectra of TbMnC>3 with the electric 
field of light polarized along each of the three crystallographic axes. We analyzed the effect, on the 
phonon spectra, of the different phase transitions occurring in this material. We show that the anti- 
ferromagnetic transition at Tn renormalizes the phonon parameters along the three directions. Our 
data indicate that the electromagnon, observed along the a direction, has an important contribution 
to the building of the dielectric constant. Only one phonon, observed along the c-axis, has anomalies 
at the ferroelectric transition. This phonon is built mostly from Mn vibrations, suggesting that Mn 
displacements are closely related to the formation of the ferroelectric order. 

PACS numbers: 75.85.+t, 63.20.-c, 78.30.-j 
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I. INTRODUCTION 

The growing interest in multiferroics was triggered by 
their potential industrial applications^ The current re- 
search focus is the puzzling mechanism driving the cou- 
pling between the magnetic and the ferroelectric (FE) 
orders, simultaneously present in these materials £~— In- 
deed, in the so-called improper multiferroics, the ferro- 
electric polarization is a consequence of a particular mag- 
netic order which breaks the inversion symmetry of the 
crystal^ In this class of materials, the magnetic origin 
of the ferroelectricity is responsible for a strong coupling 
between both orders which, for example, makes it pos- 
sible to rotate the direction of the spontaneous electric 
polarization by an applied magnetic fields 

TbMnOs belongs to this class of improper multifer- 
roics. It is paramagnetic and paraelectric at high tem- 
peratures. Below T/v ~ 41 K, the Mn spins order anti- 
ferromagnetically almost parallel to the 6-axis with an 
incommensurate wavevector k ps (0, 0.28, 1) (Ref. [y). At 
Tq ~ 28 K the Mn spins acquire a c-axis component 
which transforms the sinusoid into a cycloid magnetic 
order with approximately the same wavevector. This 
cycloid of spins breaks the inversion symmetry of the 
crystal^ and induces a spontaneous electric polarization 
along the c-axisi£ In addition to the fact that the cycloid 
magnetic and ferroelectric transitions occur at the same 
temperature, the magnetic origin of the ferroelectric or- 
der was evidenced by the fact that the electric polar- 
ization is tied to the cycloid planed and chiralityJ^ The 
electric polarization is given by P = e^+i x (S* x Sj+i) 
where S^ and Sj+i are the spins on sites i and i + 1 and 
e^i+i is the unit vector directed from site i to site i + 1. 
This magnetic spiral induced ferroelectricity was theoret- 
ically explained by Katsura, Nagaosa and Balatsky 2 - as 
an inverse Dzyaloshinsky-Moriya interaction (IDM) be- 
tween noncollinear spins creating an electronic polariza- 
tion. Inelastic neutron scattering could not detect any 



c-axis phonon anomaly related to ferroelectricity, sug- 
gesting an electronic origin for the ferroelectricityj^i An 
alternative scenario is to consider that the spiral mag- 
netic order induces ionic displacements producing the 
ferroelectric polarization. This mechanism was proposed 
by Malashevich and Vandcrbilt— whose ab-initio calcu- 
lations showed that in the ac-spiral state of TbMnO,3 
the ferroelectric polarization was mainly due to ionic dis- 
placements in the lattice. 

In TbMnOs, there is also a strong dynamic ma- 
gnetoelectric coupling which results in the existence of 
electric dipole active magnetic excitations called electro- 
magnonspi^ This electromagnon was measured by inelas- 
tic neutron scattering ) 14 i 15 Raman^ infrared (IR)^ and 
terahertz^ spectroscopies and was found to be excited by 
an electric field parallel to the a-axis regardless of the cy- 
cloid plane orientation^^ A theory based on Heisenberg 
exchange was proposed by Valdes-Aguilar et al.— to ex- 
plain this selection rule for the low frequency excitations. 
These observations suggest that the electromagnon is tied 
to the lattice rather than to the spin cycloid. Moreover, 
Takahashi et al.— evidenced a coupling between the elec- 
tromagnon and the lower energy phonon in the a-axis 
direction. For these reasons, we can expect to get use- 
ful information about the magneto-electric coupling in 
TbMnOs through the study of its lattice dynamics. 

Previous IR measurements on TbMnOs, by Schmidt 
et al.ji^ were performed along the a and c axis, but not 
the b direction. Although they have some uncertainties 
in phonon assignments due to polarization leakage, they 
found that the a and c phonon frequencies change at Tjv- 
In order to investigate more thoroughly the phonon be- 
havior of TbMnOs around the antiferromagnetic (AFM) 
and FE transitions we have undertaken a detailed tem- 
perature dependent infrared study of the phonon spec- 
tra of TbMnOs with the electric field of light polarized 
along each of the three orthorhombic directions. We were 
able to unambiguously identify the infrared active normal 
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modes predicted by group theory. The analysis of the 
thermal evolution of the mode parameters shows that the 
AFM transition affects the phonon spectra equally in all 
directions. In other words, the phonon renormalization 
at Tn along the directions (a and c) that dominate the 
dielectric properties are of the same size as the changes 
observed along the magnetic order b axis. We also find 
that along the polarization axis (c direction) additional 
changes appear at the ferroelectric phase transition. This 
is the first evidence in TbMnC>3 of a phonon anomaly re- 
lated to the ferroelectric order, adding to the debate on 
whether or not the lattice plays any role in the magneto- 
electric coupling. 



II. EXPERIMENTAL 

The samples used in this study were grown by the float- 
ing zone technique. They were aligned using Laue X-ray 
back reflection and showed dielectric anomalies at 42 K 
and 27 K corresponding to Tn and Tc, respectively^ 

We measured the near normal (10°) incidence reflec- 
tivity of two TbMnC>3 samples with their largest face 
cut, respectively, perpendicular to the 6-axis (ac face) 
and perpendicular to the a-axis (6c face). These samples 
allowed us to measure the reflectivity with the electric 
field of light polarized along each of the 3 orthorhombic 
directions at 26 different temperatures between 5 K and 
300 K in the spectral range [40-6000] cm -1 . In order to 
obtain the absolute reflectivity of the sample, we used an 
in situ gold overfilling technique.— With this technique, 
we can achieve an absolute accuracy in the reflectivity 
better than 1%, and the relative error between differ- 
ent temperatures is of the order of 0.1%. The data were 
collected in an Brukcr IFS66/s Fourier transform spec- 
trometer, the sample being attached to the cold finger of 
an ARS Hclitran cryostat. 



III. RESULTS 

TbMnC>3 crystallizes in an orthorhombically distorted 
perovskite structure with space group Pbnm {D\^). The 
irreducible representations decomposition for the normal 
modes in this system is 7A g ®8A u ®5B lg ®10B lu ®7B 2g ® 
8B 2u © 5B 3s © 10S 3u among which 9B iu , 9B 3u and !B 2u 
arc IR active for the electric field of light E parallel to the 
a, b and c axes, respectively. The 7A g , 5Bi g , 7B 2g and 
5B 3g are Raman active, all the remaining 8A U modes are 
silent. 

The solid lines on Fig. [1] are the reflectivity spectra for 
the three polarizations measured at 5 K and 300 K. The 
top panel shows the reflectivity measured with E || a. 
We find the 9Bi u modes predicted by group theory, 
which transverse optical (TO) frequencies are indicated 
by the arrows. Some features in these spectra deserve 
further discussion. Indicated by a solid down-triangle, at 
635 cm -1 there is a dip in the higher frequency phonon 
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FIG. 1. (color online) Reflectivity of TbMnOa measured with 
the electric field of light (E) polarized along each of the three 
crystallographic axis. Solid lines are the measured reflectivity 
at 5 K (blue) and 300 K (red), and the dotted lines are fits to 
the data with Lorentz oscillators. The arrows represent the 
positions of the TO phonon branches. As discussed in the 
text, the star in the top panel indicates the electromagnon 
response and the open up-triangle a possible Tb cristal field 
transition. In the top and middle panels, the solid down- 
triangles indicate LO phonon leakage due to a non normal 
incidence measurement. In the bottom panel, the solid circle 
indicates a phonon leakage due to a slight in-plane misaligne- 
ment. 



peak. This is not an additional mode of B\ u symmetry It 
is rather a signature of the highest B 3u (E |j b) longitudi- 
nal optical (LO) branch due to the fact that the incident 
light is not exactly normal to the sample's surface. This 
measurement was done with a p polarized light, which 
leads to the excitation of longitudinal modes pertain- 
ing to the crystallographic direction normal to the sur- 
face (the &-axis in the present case)<22 The peak around 
130 cm -1 , indicated by the open up-triangle, appears 



at temperatures below 150 K and was already observed 
in reflection^ and transmission^ measurements. It was 
tentatively assigned to the excitation of a double zone- 
edge magnonii, a one magnon+one phonon modei^ or a 
transition between two crystal field split levels of the Tb 
ions. At lower wavenumbers, there is a considerable evo- 
lution of the spectra below 40 K with the apparition of an 
elcctromagnon peak at 60 cm -1 , indicated by the star. 
A detailed study of the electromagnon is the subject of 
a subsequent paper. 

The middle panel of Fig. Q] shows the reflectivity mea- 
sured with E || b. In this direction, group theory predicts 
9 B3„ infrared active modes, which TO frequencies are 
indicated by the arrows. The dip at 273 cm -1 in the 
second phonon peak marked by the solid down-triangle 
is a signature, once more due to the non normal inci- 
dence measurement, of the sharp LO branch of the third 
phonon in the a direction {B\ u symmetry). 

On the bottom panel of Fig. [TJ we show the c-axis re- 
flectivity spectra at low and high temperatures. From a 
group theory analysis, we expect 7 B2 U IR active modes, 
all of them arc visible on the reflection spectrum as 
marked by arrows. The structure indicated by the solid 
circle is a polarization leakage from the a direction due 
to an in-plane misorientation of the sample, which we 
estimate to be less than 2%. The increase of the reflec- 
tivity below 75 cm -1 is an additional reflection on the 
back surface of the sample which becomes transparent in 
this region. 

In order to quantitavcly analyze the temperature evo- 
lution of the spectra, we modeled the data using one 
Lorentz oscillator for each phonon mode. In this model, 
the dielectric funtion is given by : 



e(w) =£oo +^2 



Ae k n% 



to soften due to the dilatation of the unit cell when the 
temperature is raised. This table also shows the oscilla- 
tor strengths of each phonon as well as the contribution 
of phonons and higher energy excitations to the static 
dielectric constant. The longitudinal frequencies shown 
are calculated as the zeros of the dielectric function with 
all the 7's set to zero. 
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FIG. 2. Thermal evolution of the phonon eigenfrequencies and 
oscillator strengths for E || a (Bin) between 5 K and 300 K. 
The phase transition temperatures Tn and Tc are indicated 
by the dotted lines. 



where Eqo is the contribution from electronic transitions 
to the dielectric function, and each phonon is described 
by a resonance frequency &>TO k > an oscillator strength 
Ae^, and damping 7^. These parameters were fitted so 
that the reflectivity at normal incidence, given by R — 
|1 — n/F| 2 /|l + y/e\ 2 , matches the experimental data. Fits 
are shown as dotted lines on Fig. [TJ 

The Lorentz model fits were complemented by 
Kramcrs-Kronig analysis. This allowed us to ascertain 
that very small changes observed in the fitting parame- 
ters can be tracked directly to the experimental dielectric 
function. 



IV. DISCUSSION 

The parameters of all of the IR-active phonon modes 
at low (5 K) and high (300 K) temperatures are shown 
in Table HI Except for the two lower frequency B 2u c-axis 
phonons, all eigenfrequencies are higher at low tempera- 
ture. This is the usual behavior for phonons which tend 



In order to analyze the details of the effects of the phase 
transitions on the phonon spectrum in each direction, we 
followed the evolution of the frequencies and strengths of 
each phonon with the temperature. We will not discuss 
7 values, as they not have any anomaly at Tc or Tn, 
within our fitting accuracy. 

Figure [5] shows the thermal evolution of the fitting pa- 
rameters for E || a. In this direction we observe similar 
changes as those seen by Schmidt et al.— There is a 
clear renormalization of the phonon parameters at Tn 
and within experimental accuracy, we see no changes at 
Tc- The lowest frequency phonon is the most affected 
one as its cigcnfrequcncy is increased by 3% below Tn- 
Even more striking, its oscillator strength is reduced by 
almost 50% between Tn and 5 K. These changes are 
partly due to static spin-phonon coupling, like those oc- 
curring at Tn in MnF2r^ but also to the coupling of 
the phonon spectrum with the electromagnon developing 
in the AFM phases J^ It is also worth mentioning that 
Laverdierc et al.— did not detect any spin-phonon cou- 
pling by Raman spectroscopy in TbMn03. 
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FIG. 3. Thermal evolution of the phonon eigenfrequencies and 
oscillator strengths for E || b (i?3 U ) between 5 K and 300 K. 
The phase transition temperatures Tn and Tc are indicated 
by the dotted lines. A phonon around 500 cm -1 is not shown 
as, due to its very small oscillator strength, it could not be 
fitted accurately over the whole temperature range. 
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FIG. 4. Thermal evolution of the phonon eigenfrequencies 
and oscillator strengths for Eu g ht || c (B2u) between 5 K 
and 300 K. The phase transition temperatures Tat and Tc 
are indicated by the dotted lines. Insets are a zoom at low 
temperatures for the i?2u(2) mode, which shows the largest 
change at Tc- 



The phonon parameters along the b direction, shown 
on Fig. [3] are slightly less influenced by the AFM order 
but they still exhibit clear changes at Tn- The eigenfre- 
quencies continuously increase from room temperature 
down to 100 K. Below this temperature the eigcnfrc- 
quency of each phonon mode is either further increased 
or decreased down to 5 K, surprisingly showing very lit- 
tle or no saturation. The physical response of TbMnOa 
along the b axis is dominated by the magnetic ordering 
rather than the system dielectric properties. Hence, the 
spectra renormalization observed along this direction is 
mostly due to spin-phonon coupling.^ 3 

Figure [4] shows the parameters for the c direction 
(static electric polarization axis) phonons. Although the 
most significant dielectric signatures (polarization and 
electromagnons) of TbMnOa show up in the c and a di- 
rections, the phonon renormalization along these axes are 
of the same order of magnitude as those along the mag- 
netic order b axis. Therefore, we can conclude that these 
changes are mostly due to the spin-phonon coupling. 

However our data also shows that the ferroelectric or- 
der has a clear signature along the c direction. This is 
particularly visible on the 185 cm -1 phonon mode soft- 
ening when approaching Tc from high and low temper- 
atures. As the spontaneous polarization develops along 
this axis, it is worth taking a closer look at the B2 U nor- 
mal modes. 



Figure O shows the atomic displacements for the four 
lowest frequency B 2u modes. These are the four modes 
that involve displacements of Mn atoms for this polariza- 
tion. Modes 3 and 4 are the ones that show the strongest 
changes in their oscillator strengths below Tn- Indeed, 
the atomic displacements for these modes are the ones 
that change the Mn-0 bonds the most, as the c compo- 
nent of Mn and O movements are in opposite directions. 
As a consequence, they contribute more to the renormal- 
ization of the AFM exchange (J), in a process similar to 
the one observed in MnF2»2 3 - Mode 52m (1) involves Mn 
and O displacements in the same direction. As a result, 
this phonon does not contribute to a large renormaliza- 
tion of either J or the unit cell charge distribution. This 
is in agreement with the small parameters changes seen 
for this phonon at Xjy and Tq. 

Mode i?2«(2) shows a most interesting behavior. It has 
a clear anomaly in both oscillator strength and frequency 
at Tc (see insets in Fig. 0]), indicating that Mn displace- 
ments are closely linked to the formation of the FE or- 
der. In the framework of the IDM mechanism involving 
ionic displacements, we expect a distortion of the MnC>6 
octaedra with the appearance of the FE polarization^ 
However some questions remain on whether the impor- 
tant displacements related to the FE order are in the O 
(Ref. i) or Mn (Ref. [H) atoms. As the B 2u (2) phonon 
involves mostly displacements of Mn (the oxygen atoms 
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FIG. 5. (color online) Normal modes of Bi u symmetry in- 
volving displacements of Mn atoms. 



oscillation amplitude is the smallest of the four modes 
shown) , the anomalies observed at Tq support the latter 
picture. 

We can also compare the measured static dielectric 
constant £q to the DC limit of the dielectric function ex- 
trapolated from the IR spectra. The contribution of the 
phonon and electronic excitations calculated from the fit 
to our data is given by Elf — £oo + J2 ^ £ - The static 
dielectric constant was measured by Kimura et al.— at 
10 kHz. The zero frequency infrared extrapolation gives, 



at T = 5 K (see Table U, 
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18. The values obtained by iKimura et al.l are 
e( a ) ra 24; e<» w 23 and e^ ps 29. In all cases the dc 
value is higher than the optical one. Neither the exper- 
imental errors nor sample inhomogcncities can account 
for such differences which suggests the existence of lower 
frequency excitations. Actually, a dielectric dispersion 
was observed in TbMnC>3 (and in many other RMnO,3 
compounds) along the c-axis that shows a relaxor-likc 
behavior^ Along TbMnOs c-axis, this low frequency re- 
laxation produces a drop of 10 points in the dielectric 
function when the measurement frequency increases from 
1 to 100 kHz. This drop is enough to reconcile dc and IR 
dielectric responses. Although Goto et al.— did not ob- 
serve any dispersion in the a or b dielectric constants at 
radiofrequencies, we believe that the relaxation could lie 
in the MHz or GHz regions. Dispersions in different crys- 
tallographic directions could have the same origin, once 
charge carriers usually have different activation energies 
in anisotropic systems. Even though our data alone can- 
not give a definite support the mentioned mechanism, it 
does require the existence of (sub) millimeter wave exci- 



tations. 

The thermal evolution of the low frequency dielectric 
functions along each direction can be tracked from the 
fitting parameters. Although the static values change 
along all three directions at low temperature^^ only 
the a axis data can be compared to our infrared mea- 
surements. The changes seen along b are smaller than 
our error bars and those along c are confined to a very 
narrow temperature range of 0.5 K, much smaller than 
the temperature resolution of this work. 
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FIG. 6. (color online) The symbols are the low frequency 
dielectric constant calculated from the fit results for E || a 
between 5 K and 50 K with (solid red squares) and without 
(solid black triangles) the contribution of the 60 cm -1 (HF) 
electromagnon. The dashed green line adds the estimated 
contribution of the 26 cm -1 (LF) electromagnon measured 
by Pimenov et al.— The open cirlces are the the tempera- 
tures measured by the authors. The solid blue line is the low 
frequency dielectric constant at 10kHz measured by Kimura 
et al.- 

Nevertheless, from the infrared point of view, the a 
axis dielectric constant is the most interesting as it is 
the direction where the electromagnon appears. Figure[6] 
shows the static (10 kHz) dielectric constant measured 
in the a direction by IKimura et al.l as well as the low 
frequency values extrapolated from our data. The trian- 
gles are the phonon contribution (£ p h) and the squares 
include the 60 cm -1 electromagnon (which contribution 
Aeem hf was estimated by fitting it to a Lorcntz oscil- 
lator in the reflectivity spectra). At the lowest temper- 
ature, we obtained l\eem hf ~ 3, which gives a total of 
£ph+EM HF ~ 23. Pimenov et al.— determined that a 
lower frequency electromagnon exists at 26 cm , below 
our measured spectral range. We calculated the oscillator 
stren gth (As km,*.) fo r this low frequency electromagnon 
from I Pimenov et al.l data and added it to our measure- 
ments. The result is show as the dashed line in Fig. [51 
where the open circles are the temperatures measured 



TABLE I. Lorentz fit parameters at 5 K and 300 K. Units for Qto, &lo and 7 are cm -1 . Fitted values for e^ are 5.73 along 
the a-axis 5.07 along the fe-axis and 5.37 along the c-axis. B\ u , B211 and Bs u modes are active for the electric field of light 
parallel to, respectively, a, c and b axis. The longitudinal frequencies (flho) are not input parameters to the Lorentz model. 
They were calculated by finding the zeros of the real part of the dielectric function in the absence of damping (7 = for all 
phonons) . 
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in Ref. |l3|. At the lowest temperature the contribution 
from both electromagnons solve the difference between 
static and (phonononic) dynamic dielectric constant val- 
ues. Therefore, Fig. [6] shows that the dielectric response 
is mostly due to the phonons, and the magnitude of the Eq 
jump below TV is the same as the magnitude of the elec- 
tromagnons contributions. Hence, this jump can be un- 
doubtedly assigned to the appearance of electromagnons. 
A remarkable difference occur, however, in the thermal 
evolution of the static and infrared constants between 6 
K and 50 K. Indeed, the dielectric constant at 10 kHz is 
always larger than the infrared extrapolated values, in- 
cluding the electromagnon contribution. This result is 
likely related to the observation of a dielectric dispersion 
in the c direction at radiofrcqucncies mentioned above. 
This effect, in TbMnC>3 and other related manganites, 
was tentatively attributed to the dynamics of localized 
carriers coupled to the lattice, which would affect the 
magnetic ordering^ Moreover, since the dielectric dis- 
persion exists even above T/v, one could think that fluc- 



tuations of the spin modulation wave could generate local 
regions with fluctuating net polarization, which average 
vanishes above T/v . This would imply a phase coexistence 
in a large temperature interval, as observed by Barath 
et al.— under applied magnetic field. Nonetheless, a def- 
inite explanation for the observed effect deserves further 
investigations. 



V. CONCLUSIONS 

We presented a detailed analysis of the phonon spec- 
trum of TbMnOa and its evolution with temperature. 
We showed that all the modes predicted by group theory 
arc clearly identifiable in the reflectivity spectrum of each 
orthorhombic direction. The analysis of the thermal evo- 
lution of the phonon parameters shows a clear coupling 
between the phonons in the three orthorhombic axes to 
the magnetic order. The phonon spectra renormalization 
at T/v is of the same magnitude for all crystallographic 



directions. An important result of this study is that fer- 
roelectricity produces an additional renormalization of 
the £?2 U (2) phonon along the c direction. This phonon is 
composed of Mn vibrations only, indicating that this is 
the important atomic displacement induced by the IDM 
interaction. We evidenced that in both a and b directions, 
the static dielectric constant is mostly due to the phonons 
with an additional jump at TV due to the electromagnon 
for E || a. In both cases, small sub-millimeter wave ex- 
citations are necessary to bring the dielectric function to 



its dc value. Along the c-axis the relaxor-likc behavior 
seen in Ref. ]2Q explains the much larger discrepancies 
between dc and infrared values. 
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